Broken Time Reversal Symmetry and Superconducting States in the Cuprates 
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Recently, Kaminski et al. have reported that time reversal symmetry is broken in the pseudogap 
phase in the high temperature superconducting material Bi2Sr2CaCu202+a (Bi-2212). Here we 
examine the role of translationally invariant broken time reversal states on d x i_ y i superconductors. 
In particular, we determine the change in the superconducting order parameter structure. We find 
that the broken time reversal pseudogap state that is consistent with the experiment of Kaminski et 
al., gives rise to a novel mixed singlet-triplet pairing d + ip phase. This d + ip state is shown to give 
rise to a helical superconducting phase. Consequences of this d + ip state on Josephson experiments 
■ are discussed. 

o 
o 

explanation of the pseudogap phase: it marks a new phase in which a symmetry is broken. In particular, the results of 
Ref. |J have reported that left circularly polarized photons give a different photocurrent from right-circularly polarized 
photons in the pseudogap phase. This, combined with a mirror plane symmetry, implies the breaking of time reversal 
| symmetry in the pseudogap phase. Varma 6 has proposed translationally invariant orbital current states that may 
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The origin of the pseudogap regime in the cuprates has been the subject of controversy. A variety of probes reveal 
a suppression of the single particle density of states in this regimei^. A natural explanation is that the pseudogap 
phase is a precursor superconducting state; a phase in which there are Cooper pairs but no superconducting phase 
coherence^. The recent experimental results of Kaminski et alX and of Alff et al& provide evidence for a very different 



account for the observed results. If it is indeed the case that the pseudogap phase does not break translational 
symmetry and breaks time reversal symmetry, then the classification of the superconducting pairing symmetry will 
differ from previous classifications. Fig.^shows the resulting phase diagram. In this paper, we determine the structure 
of the superconducting gap function in the pseudogap phase when time reversal symmetry is broken and translational 
symmetry is preserved. Note that this implies that we do not consider the d-density wave state of Chakravarty et alZ. 
The reason for this restriction is that a pseudogap phase which breaks translationally symmetry will not alter the 
O \ superconducting state as strongly as the case considered here. The layout of the paper is as follows: we first find the 
possible superconducting states in the pseudogap phase by using both corepresentation theory and phenomenological 
^C) I Ginzburg Landau arguments. This is done for all pseudogap symmetries that retain translational symmetry. Then 
. we focus on the pseudogap phase that is consistent with the experimental results of Kaminski et al. This phase is 
determined by the requirement that the pseudogap phase breaks time reversal symmetry, breaks the four-fold rotation 
symmetry of the Cu02 plane, and also breaks the mirror plane symmetry with normal along the Cu-0 diagonal. This 
leads to the two possible pseudogap order parameters that have been examined by Varma and Simon&S. A detailed 
symmetry analysis of the photoemission matrix elements rules out one of these two order parameters^. We develop 
■ a Ginzburg-Landau theory of the remaining pseudogap order parameter and the interplay of of this order parameter 
with superconductivity. We discuss some observable consequences of the superconducting gap components induced 
"p>! ■ by the pseudogap order parameter. In particular, we show that a helical superconducting phase is a consequence of 
broken time reversal symmetry in this case. Finally we discuss the c-axis Josephson current that has been observed 
in strongly underdoped Bi-2212 - Pb junctions^. We show how this can possibly be explained by the coupling of the 
I ' superconducting order parameter to domain walls in the pseudogap order. 

Here we describe the possible magnetic point groups of the pseudogap phase and their respective free energies. 
The symmetry group of the pseudogap phase can be written as Q—Gm X U(l) where Gm is the magnetic space 
group and U(l) is the gauge group (which is not broken in the pseudogap phase). The group Gm is the group 
that leaves both the charge density and the magnetization density M invariant. We will be interested in the pos- 
sible magnetic point groups that arise when time reversal symmetry is broken in Bi-2212. We will focus only on 
transitions that do not break translational invariance and thus focus on the 4/ mmm (D^h) point group which is 
defined through the elements {E, C 2x , C 2y , C 2z ,i, Vx,cr y ,a z ,Cda, C d b, &da, Odb, ±C iz ,±iC 4:Z }. Magnetic point groups 
are defined as Sm = H + 9(G — H), where H is a halving subgroup of the ordinary point group G. In magnetic 
symmetry groups, the crystal point group operation R, which belongs to G — H group, transforms the magnetization 
density M to — M, but the product of the time reversal operation 6 and the operation R leaves M invariant. In 
Table [IJ a list of magnetic point groups corresponding to pseudogap phase has been given. According to the exper- 
imental observations of Kapinski et al£, the four fold rotation about z-axis (C4 Z ) and the diagonal mirror planes 
(fda, Cdh), are no longer symmetry operations in the pseudogap state. The requirement that these two symmetries 
are no longer present reduces the number of possible magnetic point groups of the pseudogap phase to two: 4/ mmm 
{E, C 2x ,C 2y , C 2z ,i, a x , CT y ,CT z ,0Cda, QCdb, iOCda, Oadb, ±0C4 Z ,±i0Gi Z } or mmm {E, Cda, a z ,a d bJ9, Oada, 0C 2z ,9Cdb}. 
The symmetries, 4/mmm and mmm agree with the proposed orbital current patterns of Ref. Ifjai. Consequently, we 
use the same notation: we label the group 4/mmm as type I and the group mmm as type II. A detailed symmetry 
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FIG. 1: Possible phase diagram, temperature T as a function of hole doping, of the cuprate Bi-2212. Here T* is the critical 
temperature where time reversal symmetry is broken and T c is the transition temperature of superconducting phase. 



analysis of the photoemission matrix elements rules out the type I state. So we will consider the type II phase in 
more detail (note that Stanescu and Phillips have also examined the type II phase microscopicallyii) . 

For the one dimensional pseudogap order parameters (A or B irreducible representations (REPS) of the D^h point 
group in Table QJ, the Ginzburg Landau free energy is simply given as 

F pgJ [r 1 }=a 1 rf + ^ V \ (1) 



Two degenerate states 77 = ± yj—a\/j3\ minimize this free energy. 

For the two dimensional pseudogap order parameters (E u or E g REPS of D^h point group in Table HJl, the corre- 
sponding Ginzburg Landau free energy is 



F pgJ i[v x ,Vy] = + %) + -(Vi + %) + irf x rfy (2) 
Minimization of the free energy with respect to r\ x and r\ v gives the following sets of degenerate states: 



fe, % ) = (0,0);^{(0,±l),(±l,0)};^^{(±l,±l),(±l, T l)}. (3) 



Of these states, ?7o(±l, ±1) and ?7o(±l, Tl), where 7/0 = y—aj (f3 + 7), minimize the free energy if 7 < (3. The other- 
four degenerate states 770 (0, ±1), r/o(±l, 0), where rjo = \J —a/{3, minimize the free energy when 7 > j3. Simon and 
Varma have depicted the (1,1) state to represent the type II current pattern^, but there are three more states which 
are degenerate with this state. In these states, the tetragonal symmetry is lost as Ca z {\. 1) = (1, —1). Each domain 
has symmetry mmm with the m oriented along the different diagonals. The <Jd a (db) symmetry will be lost if there 
are multiple domains. Note that the observed incommensurate modulation in Bi-2212 implies the existence of an e xy 
strain 9 . This gives rise to an additional e xy r] x r]y invariant in the pesudogap free energy. This prefers the (±1, ±1) and 
(±1, =Fl) states and also breaks the degeneracy of these two states so that they are now two-fold degenerate. 

To classify the possible superconducting states in the presence of broken time reversal symmetry requires the use 
of corepresentation theory. The usual representation theory must be extended because the time reversal operator is 
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antilinear and antiunitary. The magnetic point group can be written as Sm = H + A x H , where A is antiunitary 
operator such that all elements of the coset AH are antiunitary. The corepresentations (COREPS) DT of Sm can 
be found from REPS T of the corresponding normal group H in one of three ways [labelled (a), (b), or (c)]. This 
approach is explained in Ref. [12| and we use their notation. Recently, similar considerations have appeared in the 
classification of superconducting states in ferromagnetsi&ilii^. The superconducting gap function is defined as^Siii 
A(k) = i(ijj(k) + id(k).&)& 2 . In Table H]we gi ye representative gap functions ip{k) and d(k) for various pseudogap 
symmetries. 

As an example, consider the 4/mmm magnetic point group. The pseudogap order parameter in this case corresponds 
to a one dimensional REP of D^. For this magnetic point group, H = D 2 h and A = 6Ci z . For a d x 2_ y 2 pairing 
symmetry when time reversal symmetry is not broken, we are interested in the Ai g REP of D 2 h- The resulting pairing 
state is a real linear combination of ^(k) = k x + k y and ip(k) = i{k x — ky) which can be denoted as a d + is pairing 
state 
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Phase(REP) 
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Basis functions 
(COREP) 


Aig(ji/mmm) 


Big 


{kl -kl)\l + ia{kl +kl)] 


A2 a (4:/mmm) 


B 2 g 


{kx ^"y) 'i>Q>kxky 


Big(4/mmm) 


Alg 


(k x — ky) + ia(k x + k y ) 


B2 g (4/ mmm) 


A 2 g 


[kx ^y ) ^^kx ky {k x ^y / 


Eg (mmm) 


Eg 


(kl - kl) + ia(rnky - r/ 2 k x )k z 


Ai u (4:/mmm) 


Bin 


(k x — ky) + ia(aik x — <J2k y ) 


j42»(4/™im) 


£>2u 


(k'i — ky) + ia(diky + (J2k x ) 


Bi„(4/mmm) 


Am 


(k x — ky) + ia(aik x + a2k y ) 


£>2u(4/mmm) 


A 2u 


(k'i — ky) + ia(diky — (J2k x ) 


E u (mmm) 


E u 


(k x — ky) + ia(f]ik y — rfekx)<f3 



TABLE I: Superconducting COREP basis functions for different magnetic point group symmetries of the pseudogap phase. 
Induced REP refers to the superconducting component induced by the pseudogap order. The magnetic point groups for two 
dimensional pseudogap REPS depend upon the form of rji, 772, here we take r/i = ±772. 

The magnetic point group mmm deserves further consideration since it is consistent with the experimental ob- 
servations of Kaminski et a&&. For this group, H = C 2v and A = i6. Due to the broken parity symmetry in the 
pseudogap phase, the pairing gap function is a mixture of spin-singlet [^(k)] and spin-triplet components [d(k)]. For 
a d x 2_ y 2 pairing symmetry when time reversal symmetry is not broken, we are interested in the Bi REP of C 2v 
(given in Table ITTf) . The corresponding COREP is a real linear combination of the spin-singlet ip(k) = (ky — k x ) and 
the spin-triplet d(k) = i(k x — k y )z gap functions. We label this state the d + ip phase. In Table II, we have also 
provided representative basis functions for the other COREPS corresponding to the different REPS of C 2v - From the 
form of the gap functions found here we can deduce whether there are possibly nodes in the gap. We find that the 
superconducting gap will vanish along the k x = k y line for in this case (the gap does not vanish along the k x = —k y 
line). Note that similar considerations appear in a recent symmetry analysis of a d-wave superconductor in a uniform 
current*^. 



Here we determine the superconducting gap structures in the pseudogap phase using Ginzburg Landau theory. 
This approach provides the same results as those found above using the less familiar corepresentation theory. Using 
corepresentation theory, it has been shown that the existence of d-wave order parameter tpd and pseudogap order 
parameter rj ensures the existence of an induced superconducing order parameter ip. In Table [IJ the different possible 
combinations of order parameters of the pseudogap order rj and of induced superconducting order ip have been 
listed. These same states can be found by examining Ginzburg Landau theory. In particular, the invariance of 
the free energy with respect to time reversal symmetry requires any free energy invariant corresponding to any one 
dimensional pseudogap REP to have the form ierj^^ip — where e is real coupling coefficient. Therefore, the 

superconducting Ginzburg Landau free energy corresponding to the one dimensional pseudogap REPS of Table [I] (up 
to second order is in ip) is : 



Fi(r], tpa, ip) = F pg j + a d \i; d \ 2 + [3 d \ip d \ 4 + a^ 2 + ientyfrl> - ^P*tpd) 



(4) 
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TABLE II: The COREPS of mmm group with H = C2v In the Table, a and b are arbitrary constants, kd a = k x + fey, 
kdb k x ky, and cri,(T2 and <T3 are Pauli matrices. The COREPS are all type (a). 

The relation between the order parameters can be obtained by minimizing F(j], ipd, ip) with respect to ip* . The relation 
is ip = Herjipd/a, where ± sign is due to the degeneracy of ip and — ip states. From the relation it is clear that for 
the non-zero r\ and ipa order parameters, ip must be non-zero, which ensures that the symmetry of superconducting 
state in pseudogap phase is d + iip. 

For a type II pseudogap phase, the order parameter belongs to the two dimensional E u REP of the D^h point 
group. Since the product of the d-wave order parameter ipd and pseudogap order parameter (r] x ,r]y) transforms as a 
E u REP, the induced superconducting order parameter also belongs to the E u REP. The induced superconducting 
order parameter can be written as -0 = (p x ,Py)- To construct a non-trivial invariant, we decompose the product of 
representations E u E u B\ g . The relevant invariant is ie[ip^(r) x p x — VyPy) ~ % l J d{nxP* x — VyPy)]i where e is a real 
positive coefficient. The Ginzburg Landau free energy [up to second order in (p x ,P v )] is 

Fn(tpd,V,P) = Fpgji + a d \ipd\ 2 + f3d\ipd\ 4 + a P (\p x \ 2 + \p y \ 2 ) + ie[ip^(rj x p x ~ f] v p y ) - tp d (VxP* x ~ VyPy)]- (5) 

Minimizing the free energy with respect to p* and p* gives p x = ih\> ] d,r] x / 1 a v and p y = —ieipd'Hy/oLp. Thus we have 
(Px,P y ) = ie(r)x, -Vy)^d/a. p = ier) (l, -l)ipd/a p where (r) x ,i] y ) = 770(1, 1) has been used. 

The appearance of the induced superconducting gap functions in the pseudogap phase leads to a variety of observable 
consequences. In the pseudogap phase having d + ip state, the breaking of parity symmetry gives rise to Lifshits 
invariants in the free energy which give rise to a spatially varying (helical) superconducting phased. This behavior 
can be readily explained in terms of the Ginzburg Landau free energy. In particular, the relevant Lifshitz invariant is 

F L = F n + ie[r d {D xPx - D v p v ) - MKpI ~ D* y p y )]. (6) 

where D = (D x ,D y ), Dj = —iVj — 2eAj/hc, and A is the vector potential. The helical superconducting phase 
can be found by setting A = and considering the spatial variation of order parameters as i\) d = V'doe lq ' r and 
iPxtPy) — {p x OiPyo) elcl ' r ■ Minimizing with respect to q x and q y gives q x = q y = — (where k is the defined through 
the gradient term kID^I 2 in the free energy). Note that gauge invariance and minimization of the free energy with 
respect to q implies that the current in helical phase is zero. The helical structure of the order parameter can be 
verified by Josephson junction experiments. We refer to Ref. |2(J for details. The existence of gap nodes found in the 
last section assumed a uniform (non-helical) order parameter. If l/q x » £0, where £0 is the coherence length, then 
the nodes will presumably still provide a reasonable description of the low energy excitations of the superconductor. 

It is interesting to note that a Josephson current has been observed through c-axis Josephson junctions between 
underdoped Bi-2212 and Pb. MoBle et al. have demonstrated this through a series of very careful experiments^. They 
have shown that the junctions are homogeneous and the coupling is a conventional lowest-order Josephson coupling. 
This experiment is difficult to explain with a pure d-wave order parameter in Bi-2212. Rae has pointed out that this 
current may exist if the Pb superconductor contains a d-w&ve contributions!. This is possible if the junction is made 
from a low-symmetry orientation of lead (a [110] face, for example). However, as pointed out by Rae, this explanation 
requires that there is a systematic bias towards [110] Pb faces in the junctions which remains to be verified. Given 
the possibility of a d + ip phase, it is natural to ask if this can be related to the observed Josephson current. Here we 
show that while a d + ip state as described above does not have a Josephson current; a domain wall in the pseudogap 
order parameter will. This can be understood by considering surface energy at a c-axis junction between a d + ip 
superconductor and a conventional s-wave superconductor; 



F sur = / d 2 S[ip*(r] x p x + rjyPy) + c.c] 



(7) 
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If we consider (r) x , rj y ) = 770(1, 1) then, as shown above, p x = —p y and F sur will be zero. But consider a domain wall 
of the type, (ri x ,r) y ) = r) [1, tanh(x/£)] (here rj x remain constant but rj y varies from -770 to 770 across the domain wall), 
the current density in this case is given as: 

j = io[l - tanh 2 (x/£)] sin(0 d - S ), (8) 

where jo = 4gt7q e | -0^ 1 1 | . A Fraunhofer pattern can be obtained by applying the field along the normal to the domain 
walls. However, as the angle between the field and domain walls decreases, the current pattern will deviate from the 
usual Fraunhofer pattern. In particular, for the field along the domain wall, only the central peak of the original 
Fraunhofer pattern remains. Note that these considerations require the junction size to be much smaller than the 
period of the helical order, so that the p-wave order parameter is approximately uniform. Furthermore, in Bi-2212, 
the incommensurate modulation may alter this analysis. If there is a strong interaction between this modulation and 
the pseudogap order parameter, then any domain walls in the pseudogap order parameter will be tied to domain 
boundaries of incommensurate modulation. 

In conclusion, we have determined the superconducting gap structure in translationally invariant pseudogap phases 
that break time reversal symmetry. It has been shown that a d + ip state is the superconducting ordered state for the 
pseudogap state that agrees with experimental results in Bi-2212. The induced ip component removes two of the four 
nodes associated with a d x i_ y i order parameter. The consequences of the induced ip order parameters on Josephson 
experiments has been explored. The induced ip phase can explain the observed Josephson current through a c-axis 
junction between underdoped Bi-2212 and Pb only if there are the domain walls in the pseudogap order parameter. 
It also has been shown that the d + ip state will give rise to a helical superconducting phase. 

We acknowledge the Donors of the American Chemical Society Petroleum Research Fund for support of this work 
and we thank V.V. Kabanov for useful correspondence. 
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